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We study the superfluid flow in a quasi-one-dimensional Fermi gas with spatially modulated inter¬ 
actions induced by an optical Feshbach resonance. Due to the competition between the periodicity 
of the modulated interaction and the nonlinearity of the background interaction, an interesting swal¬ 
lowtail structure emerges in the energy spectrum under appropriate parameters. As the interaction 
strengths are tuned, the swallowtail structure may disappear, giving rise to various different states 
on a rich phase diagram. We investigate the spatial distribution of the order parameter and the 
particle density under various parameters, which are useful for the experimental detection of the 
interesting phases in these systems. 

PACS numbers: 67.85.Lm, 03.75.Ss, 05.30.Fk 


I. INTRODUCTION 

One of the central topics in the research of cold atoms 
is the study of peculiar properties of Bose and Fermi su¬ 
perfluids. Along this line, much effort has been devoted 
to unveil various exotic features of superfluidity, e.g., col¬ 
lective excitations [1-3], vortex lattice structure [4, 5], 
the propagation of sound [6-8], and the generation of 
solitons [9]. Among them, an interesting phenomenon 
is the presence of loop structures in the energy disper¬ 
sion of a superfluid flowing within a periodic potential. 
As first discussed in the context of Bose-Einstein con¬ 
densate (BEC) [10-12], this loop, or the so-called swal¬ 
lowtail structure, is a direct consequence of the competi¬ 
tion between the nonlinear interaction and the periodic 
potential. While the periodic potential leads to a spa¬ 
tially modulated atom number density and a band struc¬ 
ture of Bloch states, the repulsive interaction between 
atoms favors a uniform distribution and tends to retain 
the quadratic dispersion. As a result, the swallowtail 
structure appears at the boundaries of the Brillouin zone 
for strong enough interactions [13]. 

For Fermi systems, the existence of swallowtail struc¬ 
tures has been discussed recently for a two-component 
Fermi gas flowing in a one dimensional (ID) optical lat¬ 
tice [14]. By tuning the interaction between the two 
spin states from the weakly interacting Bardeen-Co op er- 
Schrieffer (BCS) limit to the BEC limit via a magneti¬ 
cally tuned Feshbach resonance, a swallowtail structure 
can emerge in the energy dispersion Eg(Q ) of a super¬ 
current flowing along the lattice with a wave vector Q. 
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Although swallowtail structures in Fermi systems reflect 
different properties of superfluidity as compared to the 
Bose counterpart [14], the physics underlying the phe¬ 
nomenon is similar, i.e., the competition between the 
nonlinear interactions and the band structure induced 
by the periodical external potential. 

Besides imposing an external lattice potential, another 
route to implementing a periodical environment is via an 
optical Feshbach resonance [15-17]. In such a process, 
the coupling between the open and the closed channels 
of the Feshbach resonance is driven by a laser beam, so 
that the effective interaction between the atoms depend 
on the intensity of the laser held. If the optical Fesh¬ 
bach resonance is driven by a standing wave, the resulting 
interaction will naturally acquire a spatially modulated 
component with the same period as that of the stand¬ 
ing wave. Considering the recent experimental progress 
on optical Feshbach resonances, it is of great interest to 
investigate the effects of this spatially modulated interac¬ 
tion on the properties of a Fermi superfluid, and particu¬ 
larly, on the existence and characteristics of a swallowtail 
structure. 

In this paper, we study the properties of superfluid 
phases in a quasi-ID two-component Fermi gas with 
an interspecies interaction which is spatially modulated 
along the elongated (x) direction. The effective ID inter¬ 
action contains a uniform component, which is assumed 
to be attractive Ui < 0; and a spatially modulated por¬ 
tion I /2 cos(2fc 0 f r a:), with fc 0 f r the wave vector of the laser 
used to drive the optical Feshbach resonance. We adopt a 
mean-held approach and solve the real-space Bogoliubov- 
de Gennes (BdG) equations for pairing states with dif¬ 
ferent supercurrent wave vector Q. We And that, given 
a strong enough uniform interaction \U\\ and a weak 
enough modulated interaction U 2 , a swallowtail struc¬ 
ture can emerge in the energy dispersion Eg{Q) for Q 
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near |fc 0 f r /2|, i.e., near the boundaries of the first Bril- 
louin zone. This is consistent with the previous under¬ 
standing that the swallowtail structure is a signature of 
dominating nonlinear interactions over the periodicity of 
the system. In the opposite limit, i.e., with small \U\\ and 
large C/ 2 , the ground state of the system becomes a Fulde- 
Ferrell-Larkin-Ovchinnikov (FFLO)-like superfluid with 
a finite center-of-mass momentum Q = k 0 f r /2 [18, 19]. In 
the intermediate region between these two limiting cases, 
the ground state is either a superfluid with a critical ve¬ 
locity characterized by the critical wave vector Q c < k 0 f r , 
or a superfluid which is robust for arbitrary Q within the 
first Brillouin zone. We map out the phase diagram on 
the |Z7i I-CC 2 plane and demonstrate the spatial distribu¬ 
tions of order parameter and number distribution, which 
are helpful for the experimental detection of the different 
phases. 

The remainder of the manuscript is organized as fol¬ 
lows. In Sec. II, we introduce the model system and the 
mean-field approach we have adopted with a brief discus¬ 
sion on numerical algorithms. We discuss in detail the 
swallowtail structure in Sec. Ill, and map out the phase 
diagram in Sec. IV. In Sec. IV, we also characterize the 
FFLO-like state identified on the phase diagram. Finally, 
we summarize in Sec. V. 


II. FORMALISM 


We consider a two-component Fermi gas loaded in a 
quasi-ID confinement, where the Hamiltonian can be 
written as 


H = / dr£>t (r) 


h ^ + ^mu:l(y 2 + z 2 ) -F CT (r) 


dr 5 (r)^j(r)^j(r)^(r)>F t (r). 


( 1 ) 


Here, 'F cr (r) (dR (r)) is the annihilation (creation) op¬ 
erator for fermions with pseudo-spin (hyperfine state) 
a = (t, i) and mass m at position r. The bare inter¬ 
action g is related to the three-dimensional s-wave scat¬ 
tering length a s via the renormalization relation 1 /g = 
m/(Wh 2 a s ) — l/2ek, with the single-particle disper¬ 
sion ek = h 2 k 2 /(2m). For an optical Feslibach resonance, 
the scattering length a s ( r) can be position dependent if 
the intensity of the tuning laser acquires a spatial mod¬ 
ulation. 

In the quasi-ID limit where the transversal confine¬ 
ment frequency ui _l is much larger than other energy 
scales in the system, the motion of an atom within the 
radial plane will be frozen into the ground harmonic state 
if it is well separated from other atoms. This allows us 
to integrate out the transversal degrees of freedom and 


rewrite the Hamiltonian into an effective ID form 


H e g 


J (-^ 7 ) Vv(z) 

+ / dx-gmix^ix^l^ipiix^ix), ( 2 ) 


where ipa{x) and ip a {x) are the ID fermionic operators. 
The ID effective interaction can be related to the s- 
wave scattering length via the confinement-induced res¬ 
onance [20] 


g id 


Airh 2 a s 1 
m 7r a]j_ 



( 3 ) 


where a± = is the characteristic length of 

the transversal harmonic trapping potential, and C = 
lim s _ >oo (/ 0 ° o ds'/Vs 1 - X)s'=i I/a/s') « 1.4603. 

We then consider the case where the scattering length 
is tuned by a standing wave along the x direction, via 
an optical Feshbach resonance between the open-channel 
threshold and a molecular state in the closed channel. 
Since the molecular level is only a metastable state with 
a finite decay rate T due to spontaneous emission, the 
s-wave scattering length is typically complex 


CL S — ttbg 



Wp (I) 

H(u>e — ui c ) + ZF/2 


( 4 ) 


Here a bg is the background scattering length, the res¬ 
onance width Wq = "foI{x) is proportional to the laser 
intensity /, uii is the laser frequency, and uj c characterizes 
the detuning between the open and the closed channels 
of the Feshbach resonance. Under the practical condi¬ 
tion of Wq F, the real part of Eq. (4) dominates its 
imaginary part, and the scattering length takes an ap¬ 
proximate form in the vicinity of the resonance 


Re[a s ] 


a bg + a bg 7o-f 


- uj c ) 

K 2 {u 1 -uj c ) 2 + r 2 /4' 


( 5 ) 


This result thus leads to a spatially modulated s-wave 
scattering length a s (x ) = Re[a s ] via the position- 

dependent laser intensity I(x) = Iq cos 2 (fc 0 f r a;), where 
fcofr is the wave vector. Substituting Eq. (5) into (3), 
and focusing on the regime away from the confinement- 
induced resonance with a s < ai, we obtain an effective 
ID interaction with spatial modulations [21] 


ffiD(aO = U 1 + U 2 cos{2k oll x), 


( 6 ) 


where 


E/i = 


U 2 = 


47 rh 2 


a bg 


m TTCL , 


4 irh 2 aog 70 


7o ~ w c ) 


2 h 2 {uje -w c ) 2 +r 2 /4 

h(u>e - lo c ) 


m 2ira 2 ± h 2 (ujg — uj c ) 2 + T 2 /4 


( 7 ) 


We solve the Hamiltonian Eq. (2) via the mean-field 
approach by defining pairing order parameter A (a;) = 
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9 iu{x)(ipi(x)^(x)). The resulting Hamiltonian thus 
takes the following form 



where the chemical potential /r is assumed to be iden¬ 
tical for both spin species, the Hartree term V(x) = 
£bD(x)(4 , j(x)'I' t (x)) = gm(x)n(x) with n(x) = n t (x) = 
nj.(x) representing the spatial distribution of fermions. 
This mean-held Hamiltonian can be diagonalized with 
the Bogoliubov transformation 


Vv(x) =^2u V a(x)c va - au* CT (x)c] )(T , (9) 

v 

which leads to the BdG equations 

K(x) A(x) 

A*(x) —K{x) 

Here, c\ ] a (c Jj;tT ) is the creation (annihilation) operator of 
the quasiparticle with eigenenergy and a takes the 
value of +1 (—1) for a =f (<r =4-)- The diagonal term 
reads 

h 2 d 2 

K (x) = -^-n + v (x). ( 11 ) 

In this paper, we focus on the case of an equal spin pop¬ 
ulation, and drop the subscript of a in the solutions of 
Eq. (10) unless otherwise specified. 

The pairing order parameter and the density distribu¬ 
tion can be obtained in terms of quasiparticles via the 
following expressions 


Xrr (x) 

Xj/cr (x) 


= E, 


7](7 


Xr]c t (x) 
Xrjc t (x) 


• ( 10 ) 


A(x) = g 1D {x)'^2u ri (x)v*(x)f(E ri ), 

n( x ) = ( 12 ) 
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FIG. 1: (Color online) Energy spectra of superfluid state 
(solid lines) with supercurrent wave vector Q for (a)(c)(e) 
U^uo/Er = 0.1 and (b)(d)(f) Uiu^/Er = 0.5. The uni¬ 
form interaction \Ui\ increases from top to bottom: (a)(b) 
\Ui\n 0 /E R = 2, (c)(d) \Ui\n 0 /E R = 2.35, and (e)(f) 
\Ui\no/E r = 3, . The normal state energy En for each pa¬ 
rameter set is shown as a horizontal dotted line. All energies 
plotted in this figure are shifted relative to the superfluid state 
energy at rest Es{Q) in the corresponding cases. Notice the 
emergence of the swallowtail structures as the uniform inter¬ 
action | Hi | is strong enough. For weak uniform interaction 
| Hi |, the normal state energy is smaller than the superfluid 
state energy at Q = fc 0 f r /2, leading to a critical supercurrent 
wave vector as indicated by arrows in (a) and (b). 


over the first Brillouin zone. Denoting the total number 
of lattice sites as N L , we substitute Eq. (13) into (12), 
and solve Eq. (10) under the particle-number constraint 


Notice that the summation runs over all eigenstates of 
Eq. (10), and the function /(s) = \/(eP s — 1) denotes the 
Fermi distribution with /3 = l/ksT the inverse temper¬ 
ature. By considering a supercurrent flowing along the 
x-direction with a wave vector Q, the quasiparticle wave- 
functions can be decomposed in terms of Bloch waves 

u v (x) = e^ x J2 A ^vMx)e lkx , 

n,k 

v v (x) = e- iQx J2Bn,k,r,Mx)e ikx , (13) 

n,k 

where </> n (x) = e l2nko! * x is the wave function of the ro-th 
band, and the summation of the quasimomentum k runs 


r-n/koR 

N = N t = N± = N l / n(x)dx, (14) 

Jo 

to get the spectrum of the quasiparticles and the total 
energy of the system in the superfluid phase 

E s =J2 E vf(E v ). (15) 

v 

This energy is then compared with the normal-state en¬ 
ergy En, obtained by setting A(x) = 0, to determine the 
true ground state of the system. Numerically, we solve 
the full problem for the lowest N c bands, and employ the 
plane-wave approximation for higher bands [22]. This 




























4 



FIG. 2: Full-width of swallowtail by varying the intensity of 
the uniform interaction \Ui\. In this figure, U 2 no/E R = 0.5. 


method leads to a fairly fast convergence, and the results 
are independent of the cutoff for N c > 40. In the remain¬ 
der of this paper, we focus on the case of Ef/Er = 2.5, 
where the recoil energy Er = h 2 k 2 fT /2m is chosen as the 
energy unit. 


III. THE SWALLOWTAIL STRUCTURE 


By numerically solving the BdG equation (10) with a 
fixed filling factor n 0 = N/N L , we obtain the energy spec¬ 
tra of superfluid states as functions of the supercurrent 
wave vector Q for different interaction parameters. As 
shown in Fig. 1, the energies are relative to the value at 
rest, i.e., SEg(Q) = Eg{Q) — Eg{ 0). As the parameters 
are tuned, an interesting swallowtail structure emerges 
in a regime with a weak modulated interaction U 2 and 
a strong uniform interaction \U\\. As |Z7i| decreases, the 
swallowtail structure becomes less prominent and even¬ 
tually disappears. These results are consistent with the 
general understanding that the swallowtail structure is a 
result of dominating nonlinear interactions. 

To provide a better understanding of the swallowtail 
structures in the current system as well as the evolu¬ 
tion of the energy spectra, let us first examine the sys¬ 
tem in the weak U 2 limit, such that -C Ep (see 

Fig. l(a)(c)(e)). In this limit, the mean-field Hamiltonian 
Eq. (8) becomes spatially homogeneous to the zeroth or¬ 
der of 1 / 2 - The saddle-point solutions for the uniform 
order parameter |A(a;)| = A 0 and the chemical potential 
can be determined from the gap and the number equa¬ 
tions 


1 

V 1 


n 


_I E J_ t 8 D h(W) 

L^2E k \ 2 J 


-V 

L ^ 


1 — — tanh 

'PE k Y\ 

E k 

V 2 )\ 


(16) 

(17) 


where L is the quantization length, £& = e k — p. denotes 



k 0 fr x k 0 ‘.x 

FIG. 3: (Color online) Spatial distributions of the order pa¬ 
rameter (red) and the particle number density (blue) within 
a unit cell. Three different cases are plotted as: (a-b) SF/N 
with \U\\no/E r = 2.0; (c-d) SF with \Ui\n 0 /E R = 2.35; and 
(e-f) SW with \Ui\no/E R = 3.0. Notice in (f) that two su¬ 
perfluid states as labeled by g (ground) and e (excited) can 
be found, indicating the presence of the swallowtail structure. 
In this figure, U 2 Uo/E R = 0.5. 


the single particle dispersion shifted by the chemical po¬ 
tential, E k = + A 2 is the quasiparticle dispersion, 

and n is the number density. 

In the absence of a supercurrent (Q = 0), the total 
energy of the superfluid phase is lower than the normal 
phase for arbitrary attractive interactions (U\ < 0). The 
condensation energy Ec, which is defined as the energy 
difference between the normal state En and the super- 
fluid phase Eg., increases with \U\\. By imposing a su¬ 
percurrent with finite wave vector Q 1 the total energy of 
the superfluid acquires an addition kinetic term 


SE s (Q) 


= E S (Q)~E S ( 0) 

ft 2 Q 2 A _ _£fc_ 

2m k A E k 


tanh 


pE k 

2 


21 k^k •> 
(18) 


where Y k = ( f3/2)sech 2 {j3E k /2) is the Yoshida distribu¬ 
tion. The critical supercurrent wave vector is then de¬ 
fined as the critical Q c where the gained condensation 
energy gets compensated by the kinetic energy associated 
with the superfluid flow, i.e., Eq = SEs(Q c )■ When the 
superfluid flows faster than the critical Q Cl the superfluid 
ground state undergoes a phase transition in the ground 
state to become normal. 
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FIG. 4: (Color online) Phase diagram for the parameter re¬ 
gion Ui < 0 and U 2 > 0. Four regions discussed in texts can 
be identified on this diagram, and are separated by first-order 
(solid) or second-order (dashed) transition lines. 


The inclusion of a small interaction term U 2 will induce 
small spatially periodic modulations in both the number- 
density and the order-parameter distributions on top of 
a uniform background. This gives rise to a periodic en¬ 
ergy spectrum in the quasimomentum space. As a con¬ 
sequence, the energy of a superfluid flow with quasimo¬ 
mentum Q also acquires a periodic structure with multi¬ 
ple copies in higher Brillouin zones, with different copies 
intersecting at positions of Q = (£ + l)fc 0 f r /2 (£ is an 
integer). For a weak uniform interaction \U\\, the con¬ 
densation energy is smaller than the kinetic energy at 
the intersections, and the supercurrent critical wave vec¬ 
tor Q c is naturally smaller than k a f r /2, as indicated by 
arrows in Fig. 1(a). In this case, the system is a super¬ 
fluid with conventional dispersion for slow supercurrent 
flow (| Q\ < Q c ), and becomes a normal fluid with fast 
flow (|Q| > Q c ). We denote such a case as SF/N. On the 
other hand, when the uniform interaction | U\ | becomes 
stronger, the condensation energy can be so large that 
the normal state energy lies above SEs(k 0 [ r /2). In this 
case, the swallowtail structure can emerge, under appro¬ 
priate interaction parameters, as shown in Fig. 1(e). 

As the modulated interaction strength U 2 increases, 
the spatial modulations of number density and order pa¬ 
rameter become more evident, and provide a stronger 
periodic environment for the particles. The strong spa¬ 
tial periodicity of the system thus leads to an apparent 
band structure in the quasimomentum space with a gap 
opening at the boundaries of the Brillouin zone. In other 
words, if we consider the superfluid phase with supercur¬ 
rent wave vector Q as a macroscopic quantum state, the 
presence of an interaction term U 2 can induce an effec¬ 
tive coupling between supercurrents at different bands, 
such that a gap can be opened at the location where the 
two bands intersect. In this case, if the uniform interac¬ 
tion \Ui\ is weak such that the normal state energy lies 
within the lowest band, the system is in SF/N state as we 





FIG. 5: (Color online) (a) Energy spectrum of the superfluid 
state with supercurrent wave vector Q for modulated interac¬ 
tion strength U 2 Tlo/Er = 3.5. Notice that for small enough 
uniform interaction |t/i|, the superfluid state with Q = k 0 f T /2 
can become energetically favorable, so that the ground state 
of the system becomes an FFLO state, (b) Spatial distri¬ 
butions of the superfluid order parameter A(x) (red, solid) 
and the particle number density n[x) (blue, dotted) for the 
FFLO state within a unit cell. Parameter used in this figure 
is Uino/ER = 3.5. 

have discussed previously (see Fig. 1(b)). Conversely, if 
the uniform interaction is so large that the normal state 
energy lies above the lowest band, the ground state of 
the system is always a superfluid (SF) (see Fig. 1(d)). If 
\Ui\ further increases, the uniform interaction overcomes 
the effect induced by U 2 and restores the quadratic form 
of ground state energy. As a consequence, the system 
acquires the swallowtail structure for Q around fc 0 fr/2, 
as indicated in Fig. 1(f). The width of the swallowtail 
Qsw for a fixed U 2 is shown in Fig. 2, from which one 
can see clearly that the swallowtail extends with increas¬ 
ing \U\\. This observation is consistent with the under¬ 
standing that the uniform interaction favors swallowtail 
structure. 

The real-space distributions of number density and or¬ 
der parameter for the three different cases are listed in 
Fig. 3, using the same parameters as in Fig. l(b)(d)(f). 
Notice in Fig. 3(d) that the spatial distributions show 
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clear double-peak structures within a unit cell for the 
case of Q = k 0 f r /2 within the SF case. This is a clear 
signature of the higher band effect induced by the inter¬ 
action term [/ 2 at the gap-opening position. 

IV. PHASE DIAGRAM AND THE FFLO-LIKE 
STATE 

With these understandings, we map out the phase di¬ 
agram for the parameter region Ui < 0 and t/ 2 > 0. As 
shown in Fig. 4, the swallowtail structure can be found in 
a sizable parameter region (SW). Additionally, the SF/N 
and the SF regions can also be identified, consistent with 
our previous analysis. At small \U\\ and large U 2 how¬ 
ever, an FFLO-like state with a finite center-of-mass wave 
vector Q = k Q f r /2 emerges on the phase diagram. The ex¬ 
istence of this FFLO-like state can be understood along 
the same line of the previous energy spectrum analysis. 

Starting from the SF/N region, a decreasing \U±\ leads 
to a weaker non-linear interaction, or equivalently, a 
stronger modulation in space. As a consequence, the pe¬ 
riodicity effect becomes more evident, which enhances 
the effective coupling between the two bands, and leads 
to the gap opening between the two superfluid states at 
Q = k 0 fr/2. The superfluid energy E$ of the lowest band 
at Q = fcofr/2 thus decreases, as depicted in Fig. 5(a). 
As the energy Es{k 0 f r /2) becomes lower than the energy 
Es{ 0) of the superfluid state at rest, the system under¬ 
goes a first-order phase transition to become a super¬ 
fluid state with a spontaneously generated finite center- 
of-mass momentum, i.e., an FFLO-like state [18, 19]. The 
phase boundary of this FFLO-like state is determined by 
the position at which Es(k 0 f I /2) and E, g(0) are degen¬ 
erate (solid line in Fig. 5(a)). In Fig. 5(b), we show 
typical spatial distributions of the number density and 
order parameter for the FFLO-like state. Notice that 
due to the strong modulated interaction t/ 2 , the order 
parameter acquires nodes within a unit cell. Besides, the 
two degenerate states with Q = ±fc 0 f r /2 spontaneously 
break inversion symmetry, as can be usually observed in 
a single component Fulde-Ferrell state [18]. As a conse¬ 


quence, while the SF region is separated from SF/N and 
SW by second order phase boundaries, the FFLO-SF and 
FFLO-SF/N transition lines are of first order because of 
the abrupt change of the ground state from Q = 0 to 
Q = k 0 { T /2 at the phase boundary. 


V. CONCLUSIONS 

We investigate the superfluid states in a quasi-ID 
Fermi gas under a spatially modulated interaction, which 
can be realized by an optically tuned Feshbacli resonance. 
By calculating the ground state energy of the system 
via a BdG approach, we obtain the energy spectrum of 
a superfluid flowing with wave vector Q, and reveal a 
swallowtail structure near the boundary of the Brillouin 
zone within the regime of strong attractive background 
interactions and weak modulated interactions. This is 
consistent with the physical picture that the swallowtail 
structure is present as the effect of a nonlinear interaction 
overcomes that induced by the periodic environment. We 
map out the phase diagram by varying the interaction, 
and identify four different regions therein, including an 
interesting FFLO-like region. For each region, we inves¬ 
tigate typical spatial distributions of the order parameter 
and particle density, which are helpful for the experimen¬ 
tal detection. 
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